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We present the results of our observations of eight magnetic Herbig Ae/Be stars obtained with the X-shooter spectrograph 
mounted on UT2 at the VLT. X-shooter provides a simultaneous, medium-resolution and high-sensitivity spectrum over 
the entire wavelength range from 300 to 2500 nm. We estimate the mass accretion rates (M acc ) of the targets from 13 
different spectral diagnostics using empiric calibrations derived previously for TTauri-type stars and brown dwarfs. We 
have estimated the mass accretion rates of our targets, which range from 2 x 10 -9 to 2 x 10 -7 Mq/jt. Furthermore, we 
have found accretion rate variability with amplitudes of 0. 10-0.40 dex taking place on time scales from one day to tens of 
days. Additional future night-to-night observations need to be carried out to investigate the character of M acc variability 
in details. Our study shows that the majority of the calibrational relations can be applied to Herbig Ae/Be stars, but several 
of them need to be re-calibrated on the basis of new spectral data for a larger number of Herbig Ae/Be stars. 
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1 Introduction 

Herbig Ae/Be stars (HAeBes) are pre-main sequence (PMS) 
objects of intermediate mass approximately from 2 to 
10 M Q (Herbig [T9601 Finkenzeller & Mundt [l98"4i l. They 
range from spectral class F2-09 and are surrounded by 
relict dust/gas accretion disks, reflecting themselves in the 
form of an IR excess as a result of thermal emission of 
circumstellar (CS) dust. For a number of stars, the disks 
can also be directly detected on images (e.g. Grady et al. 
2005 ). The CS environment of HAeBes possesses a com- 
plex structure. Gaseous streams are accreted from the equa- 
torial disk onto the star and a stellar wind zone exists at 
higher latitudes. On the other hand, the disk wind, the so- 
called magnetocentrifugal wind, carries the excess of angu- 
lar momentum away (e.g. Garcia et al. 2001 ). CS contribu- 
tion to the stellar atmospheric spectrum provides informa- 
tion on physical processes in the envelope and on its inter- 
action with the star. The character of this interaction is rea- 
sonably well studied for PMS objects with masses of about 
one solar mass and less (TTauri-type stars or TTS). These 
objects have significant (several kG) dipole-type magnetic 
fields. The star/CS interaction in TTS is described in de- 
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tail by the model of magnetospheric accretion (MA; e.g. 
Bouvier et al. 2007 and references therein). According to 
the MA-model, the field truncates the disk where the ram 
pressure of the accreted matter is balanced by pressure from 
the stellar magnetic field. At the truncation radius the ac- 
creted material is channeled into funnel flows and ballisti- 
cally falls onto the stellar surface. A situation that is less 
clear for HAeBes. Their interiors are mostly radiative (see 
e.g. Fig. 8 in Hubrig et al. [2009). However, convection is 
needed to maintain a stellar dynamo that generates a strong 
magnetic field. Nevertheless, most late-type HAeBes with 
spectral classes later than B9 (HAes), have possibly evolved 
from the so-called intermediate mass T Tauri stars (IMTTS) 
of K-G spectral types (see e.g. Palla & Stahler 1993) pos- 
sessing sub-atmospheric convective zones and, as a result, 
strong magnetic fields. It is likely that these fields do not dis- 
sipate completely during the PMS evolutionary stage. Ad- 
ditionally, it is possible that magnetic fields of HAeBes are 
generated by differential rotation of the star or originate in 
dynamos within the CS disk (Tout & Pringle 1995; Tout 
& Pringle |19961 l. Recent direct measurements of HAeBe 
magnetic fields using spectropolarimetry (Wade et al. 2007 
Hubrig et al. |2007l 125091 1201 lal and references therein) 
have allowed to detect longitudinal magnetic fields of hun- 
dreds of Gauss in a number of HAes. Stars for which so far 
no magnetic field could be detected may still have a field 
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but below the current sensitivity limits. The inner boundary 
of their accretion disks are likely to be located close to the 
stellar surface. 

The existence of magnetically controlled disk accretion 
in HAes going through polar funnels was confirmed by 
Muzerolle et al. (120041 1. Mottram et al. (120071 1. and Grady 
et al. (1201 Oi l. In any case, the magnetospheric radii of HAes 
with BwlOOG must be much smaller in comparison with 
those of TTS, where magnetic fields are of the order of 
several kG. As a result, we can expect the character of the 
star/CS interaction in HAes to be different from that of TTS . 
In our study aim at investigating the general picture of the 
star/CS interaction in HAes and at comparing it with that in 
TTS. One of the initial goals was also to study the vari- 
ations of the accretion mass rates over the magnetic pe- 
riod. However, due to weather or technical constraints the 
observations were carried out mostly for targets with sin- 
gle magnetic field measurements, i.e. for targets with un- 
known magnetic phase curves. Clearly, an analysis of the 
temporal behaviour of different mass accretion rate indica- 
tors over the magnetic/rotation periods would provide fun- 
damental clues on how magnetospheric accretion works. 
The targets for which magnetic/rotation periods have been 
determined in our previous studies include HD 97048 and 
HD 101412 with only one X-shooter observations, respec- 
tively, HD 150193 with two observations, and HD 176386 
with four observations (Hubrig et al. 1201 lb] 1201 la) . 

One of the most important parameters characterizing the 
properties of the star/CS interaction is the mass accretion 
rate Af acc . It determines strongly the disk structure and dy- 
namics (e.g., D'Alessio et al. 1999, 2001). Measurements 
of Af acc are well-developed for TTS, based on a number of 
empirical spectral diagnostics (e.g., Oudmaijer et al. 1201 II 
Rigliaco et al. 1201 II and references therein). 

Three main purposes of our work presented here are: 

1 . To determine A/ acc in a sample of HAes whose mag- 
netic fields have been detected (Hubrig et al. 120091 
1201 lbl l using calibrations previously obtained for TTS 
and brown dwarfs (BDs). Magnetic HAes appear to be 
most suitable targets for such a study because their disk 
accretion process is expected to better fit the MA sce- 
nario in comparison with non-magnetic stars. 

2. To test the applicability of different calibrations for mea- 
suring Aiacc of HAes and to select the best spectroscopic 
tracers of accretion necessary to provide accurate ac- 
cretion rates for systematic studies of large samples of 
HAes. 

3. To search for variability of Af acc in HAes with more 
than one available spectrum. This variability can be 
twofold: a) a real change in the accretion process, and 
b) a modulation of Af acc by the rotating magnetosphere 
and a cyclic screening of the stellar limb by local funnel 
streams. 



2 Mass accretion rate calibrations 

2.1 Calibrations for TTS and BDs 

At present, a number of empiric spectral indicators are avail- 
able for the determination of M acc for TTS and BDs. 

2.1.1 Model method 

The basic accretion indicator is the additional emission in 
the continuum superimposed onto the stellar spectrum and 
calculated on the basis of MA models for these objects 
(e.g., Calvet & Gullbring [1998l Gullbring et al. EOOOl Ardila 
et al. I2002I I. It is assumed that the overwhelming bulk of 
this emission originates in funnel flows inside the magneto- 
sphere and in the shock zone on the stellar surface in a com- 
pact region. This assumption is rather plausible since the 
emissivity in the continuum is proportional to N% (where 
N e is the electron density), and the magnetospheric streams 
consist of the densest CS matter. The gas density in other 
parts of the envelope (disk and wind) drops very quickly 
with distance from the star, and its contribution to the ad- 
ditional continuum emission is negligible. In the observed 
spectra, this emission is seen as veiling of the atmospheric 
lines and as an excess radiation at wavelengths below the 
Balmer jump. The theoretical models allow to construct a 
quantitative relation between A/ acc and the veiling plus UV 
excess. Yet, veiling and UV excess frequently prove not to 
be convenient for these measurements, due to the large ex- 
tinction in the regions where the targets are located. 

2.1.2 L - type calibrations 

Alternatively, empirical calibrations have been determined 
between emission lines and mass accretion. 

All emission lines in the spectrum of a PMS object are 
of multi-component origin. They are formed not only in 
the magnetospheric streams, but also in the disk and in the 
wind. The source function determining the emissivity in a 
line drops with the distance from the star much slower than 
the emissivity in the continuum (for example, see Pogodin 
1986 1989 for Balmer emission lines). However, according 
to the MA models, such parameters as the mass of the star, 
mass of the disk, A/ acc and mass loss rate Af w j nc ; are likely 
to be interdependent. It has been established that the accre- 
tion luminosity (L acc ) is clearly correlated with luminosity 
in some emission lines (Lii nc ). A number of empirical cali- 
brations have been published so far . The dependencies are 
approximated by expressions of the type: 

l0g(L a cc) = alog(Llinc) + b, (1) 



where Af* and are the stellar mass and radius, re- 
spectively, and a and b are coefficients determined sepa- 
rately for each spectral line (e.g., Fang et al. [2009. Dahm 
2008 Herczeg & Hillenbrand |2008l Gatti et al. |2008l Natta 
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Fig. 1 Relationship between ADb and M acc constructed 
on the base of data extracted from Table 3 in Donehew & 
Brittain ( 1201 11 1 and fitted by a polynomial. 

et al. 120041 Mohanty et al. l2.005t . The calibrations are con- 
structed on the basis of spectroscopic observations of the 
chosen line in a representative sample of targets, whose 
M acc values were determined beforehand using the model 
method. 

2.1.3 F - type calibrations 

Besides of that, another type of empirical calibrations ex- 
ists that is based on a direct correlation between M acc and 
the flux in the emission line (-Fline). The dependencies are 
expressed as: 

log(M acc ) = clog(Fu ne ) + d, (3) 

where c and d are coefficients (e.g. Herczeg & Hillenbrand 
2008). The F-type calibrations are constructed analogous to 
the case of L-type calibrations where the values of M acc are 
derived beforehand using the model method or from other 
types of calibration. 

2.1.4 Hal0% - indicator 

There is also a separate accretion indicator, the Hal0%. It 
is defined as the full width W of the emission Ha profile at 
the level of 10% of the maximum intensity (White & Basri 
12003b . Natta et al. d20041 > extended this indicator to very 
low-mass objects, where M acc was determined by fitting 
the observed Ha profiles with predictions of MA models 
(Muzerolle et al. 120011 1. The calibration is expressed in the 
form: 

log(M acc ) = d \og(W(Hal0%)) + di, (4) 

where c\ and d\ are coefficients and W(Hal0%) is given in 
(km/s) (Natta etal.EOQl 



2.2 Calibrations for HAeBes 

At present, calibrations constructed specifically for HAeBes 
are small in number. 

Donehew & Brittain (2011) derived M acc for a consid- 
erable sample of HAeBes using the model method, where 
the modeling was applied to the emission Balmer jump 
ADb - They calculated M acc in dependence of the emission 
Balmer jump for more than 30 HAeBes using the models 
and methods developed by Calvet & Gullbring (1998), Gull- 
bring et al. (l2"0"0"0i l. and Muzerolle et al. (l2"004i l. The model 
of Muzerolle et al. (2004) takes into account the main phys- 
ical processes predicted by the MA model such as: a bal- 
listical infall of material from the disk onto the stellar sur- 
face along accretion columns, heating of the photosphere 
due to its shock interaction with the accretion streams and 
releasing soft X-rays that are absorbed by the surrounding 
material. This material then emits optical and UV radiation 
as it thermalizes. In the model, the flux from accretion is 
calculated for three different regions: the shock region, the 
heated photosphere, and the accretion columns, to get the 
overall flux from accretion, which is a function of M* , i?* , 
and Macc- The resulting spectral energy distribution (SED) 
including both the accretion and the stellar flux is used to 
calculate the expected emission Balmer jump ADb for dif- 
ferent M acc and stars of different spectral types. The value 
of ADb was defined as the difference in Db (Balmer jump) 
for a given star from that of a standard (non-accreting) star 
of the same spectral type, where ADb is the difference in 
the magnitudes at both sides of the discontinuity. The rela- 
tionship between ADb and L acc calculated by Muzerolle et 
al. ((2004) and between L acc and Af acc given in Sect. 12.1.21 
were used by Donehew & Brittain (1201 1) to determine M acc 
of their program stars with different M* and i?* . They used 
the expression: 

AD B = 2.b\og F * +FA . (5) 

-T* 

The observational AD b for all targets have been mea- 
sured in the observed spectra normalized to the flux F\ at 
4000 A using the spectra of standard stars and the Castelli- 
Kurucz LTE models (Castelli & Kurucz 2008 ) of stellar at- 
mospheres with corresponding T ff and log g. As it turned 
out, the observational ADb demonstrated a very convinc- 
ing correlation with theoretical M acc . This result is likely to 
be related to the fact that the MJR* -ratio is practically the 
same in PMS stars from late-B to K types with a dispersion 
of ±0.1 dex. We have built the dependence between ADb 
and M acc using the data taken from Columns 2 and 3 of Ta- 
ble 3 in Donehew & Brittain (1201 11 1 and fitting them by a 
polynomial. This relation is shown in Fig. [TJof this paper 
and used in our study. 

Donehew & Brittain (1201 11 1 measured also the equiv- 
alent widths of the B1-7 emission line and constructed the 
empiric correlation between log(i acc ) (determined from 
ADb) and Lb^- They concluded that this relationship is 
different for HAes and for earlier type HAeBes (HBes) 
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and that the correlation for HAes is in satisfactory agree- 
ment with the calibration constructed earlier by Calvet et 
al. (120041 for classical TTS (CTTS) and IMTTS. Therefore, 
it appears that the B1-7 emission line is probably a reliable 
tracer of Af acc in HAes. One of the goals of our study is test- 
ing the reliability of additional available spectral indicators. 

Also Garcia Lopez et al. (120061 1 determined Af acc for 
36 HAes using a similar calibration for TTS. Eight objects 
are in common between their target list and the list of Done- 
hew & Brittain ( 1201 11 1. For four objects the agreement in the 
derived Af acc values is rather good with A log Af acc rang- 
ing from 0.01 to 0.23 dex. Four other objects demonstrate 
significant differences from 0.7 to 1.2 dex. As a whole, the 
difference between the results of these two works makes up 
—0.18 ± 0.75 dex, where the error is the standard deviation. 
Taking into account that the accuracy of an individual es- 
timate is of the order of ±0.4-0.5 dex in both papers, the 
conclusion can be made that the results obtained in these 
two works demonstrate no significant difference. 

During the time when our work was in preparation 
one additional paper was published by Mendigutia et al. 
(1201 lal l. where a similar model method was used to deter- 
mine mass accretion rates in 38 HAeBes and the L-type cal- 
ibrational relationships were constructed for several emis- 
sion lines. We compared the results of this work with those 
obtained in Donehew & Brittain (1201 11 1 and Garcia Lopez 
et al. ( 120061 1. The Mendigutia et al. ( 1201 lal l sample has 
overlap with the target list of Donehew & Brittain (1201 1) 
- 14 stars - and with that of Garcia Lopez et al. {2006i 
-13 stars. We see that Af acc of about half the targets of 
Mendigutia et al. d201 lab are estimated to be systematically 
higher in comparison with those obtained in the two other 
studies (+0.74 ± 0.44 dex for Donehew & Brittain |20TT1 
and +1.67 ± 0.83 dex for Garcia Lopez et al. 120061 Only 
upper limits of M acc have been obtained in Mendigutia et 
al. (1201 lai for the second half of their targets. The differ- 
ences in their Af acc determinations with those obtained in 
the two other works account for more than — . 74 ± . 56 dex 
compared to the work of Donehew & Brittain ( 201 11 1 and 
—0.50 ± 0.24 dex compared to the work of Garcia Lopez et 
al. (l2006h . 

Trying to understand the possible cause of these dis- 
crepancies, we compared the values of the emission Balmer 
jump ADb measured in Mendigutia et al. ( 1201 lal l by a pho- 
tometric method and in Donehew & Brittain (1201 11 . where 
a spectroscopic method was used. We have obtained the dif- 
ference 0.17 ± 0.17. Taking into account the different dates 
of observations and different methods of measurements and 
their errors (±0.07 for Mendigutia et al. 120 11 al and ±0.03 
for Donehew & Brittain 1201 11 1, we can conclude that the 
differences in these measurements are not essential. 

Thus we assume that some systematic distinctions can 
be expected in a model calculating the relationship between 
Af acc and ADb carried out in the different studies. We note 
that Af acc values determined in Mendigutia et al. (1201 lal l 
appear significantly overestimated: ten objects in their study 



Table 1 Target stars for which X-shooter data were ob- 
tained during our observing run. Spectral type and photo- 
metric data were taken from the SIMBAD database. 



Object 


Spectral 
type 


V 


R 


J 


H 


A v 


HD 97048 


AOpshe 


8.46 


8.50 


7.27 


6.67 


1.00 


HD 100546 


B9Vne 


6.80 


6.70 


6.43 


5.96 


0.10 


HD 101412 


B9.5V 


9.29 


9.30 


8.64 


8.22 


0.54 


HD 135344B 


F4-F8 


7.85 


7.83 


7.31 


6.67 


0.10 


HD 150193 


AlVe 


8.88 


8.90 


6.95 


6.21 


1.60 


HD 176386 


B9IVe 


7.30 




6.90 


6.75 


0.60 


HD 190073 


A2IVpe 


7.82 


7.80 


7.19 


6.65 


0.12 


PDS2 


F2 


10.73 




10.01 


9.68 


0.55 



present rates ranging from 10~ 4 to 10~ 6 AfQ/yr. Such large 
mass accretion rates can be expected only in protostars or in 
such unique objects as FUORs. In addition, according to the 
stellar masses and ages presented in Table 1 of Mendigutia 
et al. (1201 lal l, six out of ten targets of the spectral type B8 
- A2 are located on the HR diagram close to the main se- 
quence where the accretion rates already become smaller 
than at earlier stages of the PMS evolution. Generally, the 
masses of accretion disks around HAeBes range from 0.01 
to O.IMq (Hillenbrand et al. [19981 Henning et al. [19981 
Natta et al. 120001 1. With mass accretion rates of the order 
of 10~ 4 to 10 -6 AfQ/yr, the accretion disk would be com- 
pletely dissipated already after 10 2 — 10 5 years. This time 
interval is much smaller than the time of the whole PMS 
evolution stage of objects with masses of 2Mq, which is 
about 10 6 years according to Palla & Stahler (119931 . 

Based on these results, the question arises whether cal- 
ibrations determined for lower mass pre-main sequence 
stars, TTauri stars, can be employed to estimate A/ acc in 
HAeBes. Clearly, applicability of such calibrations for Af acc 
in HAeBes demands a special examination. It can be as- 
sumed that identical calibrations are applicable to HAeBes 
provided that the character of the disk/magnetosphere and 
the disk/wind interaction in TTS and HAeBes is similar and 
that the contribution of the CS regions to the whole CS spec- 
trum formed in TTS and HAeBes is comparable too. As was 
discussed in Sect. [T] the size of the magnetosphere and the 
value of B is expected to be much smaller in HAeBes in 
comparison with TTS. It is not clear yet whether these dif- 
ferences have an impact on the relations between Af acc and 
the emerging radiation in spectral lines. 

3 Observations and data analysis 

The observations were performed in service mode with the 
X-shooter spectrograph mounted on the 8 m UT2 of the 
VLT. X-shooter allows to obtain spectral data simultane- 
ously over the entire wavelength range from the near-UV to 
the near-IR in three arms (the UVB-arm covering the range 
300-590nm, the VIS-arm 550-1000nm, and the NIR-arm 
1.0-2.5 /j,m). The observations were performed with the 
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Table 2 Stellar parameters of the targets. The majority of the data are taken from Hubrig et al. ([2009 1. Further data is 
from; a - Catala et al. ([2007 ); b - Montesinos et al. (2009); c - Coulson & Walter ( 119951 ); d - this work. 



Object T cff log g M./Mq R,/Rq usini < B z > 

[K] km/s [G] 



HD 97048 


10000 


4.0 


2.5 


2.2 


140 


188±47 


HD 100546 


10500 


4.5 


2.5 


1.5 


65 


89±26 


HD 101412 


9500 


4.0 


2.5 


2.7 


5 


-454±42 


HD 135344B 


6250 c 


4.0 C 


1.36 


1.9 


70 


-37±12 


HD 150193 


9000 


4.0 


2.2 


2.0 


100 


-252±48 


HD 176386 


10000 


4.0 


2.7 


2.3 




-121±35 


HD 190073 


9250 


3.5 


2.85" 


3.6° 
2.1 6 


12 


104±19 


PDS2 


7000 d 


4.0 d 


2.5 


1.6 


30 d 


103±29 



Table 3 Non-accreting stars used as flux standards 



Object 


Spectral type 


MJD 


V 


HD 100604 


F2V 


55352.008 


7.69 


HD 100627 


F6IV/V 


55408.332 


8.49 


HD 100926 


A3III/IV 


55466.170 


9.75 


HD 100928 


A0IV 


55648.147 


9.52 



highest possible spectral resolution, i.e. R is —9100 in the 
UVB-arm, -17400 in the VIS-a rm, and -11000 in the 
NIR-arm (D'Odorico et al. [2006 ). The data were reduced 
using the X-shooter pipeline (Version 1.1.0) following the 
standard steps. For more details see Modigliani et al. (120 101 ). 
Due to a very high efficiency of the X-shooter spectrograph 
the signal-to-noise ratio (S/N) of 300-500 was achieved 
during exposure times ranging from 13-15 s for the bright- 
est targets to 450 s for the faintest target PDS 2. 

26 spectra of eight HAes with magnetic field detections 
were obtained during 13 nights distributed between May 
and September 2010. After each science exposure, telluric 
standards were automatically observed in Obseratory time 
at S/N— 100. They are usually main-sequence hot stars or 
solar analogs. Further, the data package delivered by ESO 
includes a number of spectrophotometric standards. The 
flux standards used in our work are presented in Table |3] 

The original observation request foresaw several obser- 
vations per target over the rotation period, but due to tech- 
nical problems with X-shooter and mediocre weather con- 
ditions, the program was only partially completed. The tar- 
gets were chosen from the sample of Herbig Ae/Be stars 
investigated previously by (Hubrig et al. 12009) . The list of 
X-shooter targets is presented in Table [TJ and their stellar 
parameters together with the detected longitudinal magnetic 
fields are summarized in Table [2] Since the star PDS 2 was 
less intensively studied in the past than other targets, we 
estimated atmospheric parameters of this object by compar- 
ison between the observed spectrum in the range AA4460- 
4500 A and the synthetic spectrum computed using the code 
SYNTH+ROTATE (Piskunov[[992]). The best model fit is 
presented in Fig. [2] with parameters: T e s= 7000 K , logg = 



1.1 

1.0 
0.9 
0.8 
0.7 
0.6 

1.0 
0.9 
0.8 
0.7 
0.6 



Model: 

T,„ = 7000°K, logg= 4.0. V sini = 30 km/s 




PDS2 (mean) 



4460 



4470 



4480 

A (A) 



4490 



Fig. 2 Comparison of the mean observed spectrum of 
PDS 2 in the range of AA4460^4500A with the synthetic 
spectrum corresponding to the best fit. Top panel: The syn- 
thetic spectrum. Bottom panel: The overplotted observed 
spectrum (solid line) and the synthetic spectrum (dotted 
line). 



4.0, and vsini — 30 km/s. The complete list of observing 
dates is given in Table[5] 

The data analysis concentrated on the emission lines and 
the emission Balmer jump that can be considered as accre- 
tion indicators. To remove telluric features, we used the tel- 
luric standards observed immediately after each observation 
at a similar zenith distance. The equivalent widths (EWs) of 
emission lines were measured after subtraction of stellar at- 
mospheric profiles from the observed ones. The synthetic at- 
mospheric profiles were calculated using the computer code 
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SYNTH+ROTATE (Piskunov|T992| and the standard LTE 
Kurucz models for the corresponding values of the stellar 
parameters (see Table |2}. The synthetic atmospheric line 
profiles were also smoothed in accordance with the spec- 
tral resolution and vsini. For the near-IR hydrogen lines 
(Br7, Pa/3, and Pa7), the excess of radiation was taken into 
account by comparing the photometric magnitudes of the 
target stars with those for unreddened stars of correspond- 
ing spectral type and a correction for A\ using the standard 
extinction law. All necessary data for this step are presented 
in TableQ] The EWs of the Na I D lines were measured as a 
sum value for both components D\ and D-2, where the wa- 
ter vapour absorption lines and the interstellar (IS) narrow 
components were removed prior to this step. The IS com- 
ponents of the Na I D lines were cleaned by simple cutting. 
This procedure leads to a small underestimation of the EW, 
but it does not account for more than 10% (~0.1 dex) in 
the value of logM acc . The emission lines of the IR Ca II 
triplet are blended by weak emission components of Pal 3, 
Pal5, and Pal6. But their contribution to the EWs of the 
Ca II lines can be taken into account by polynomial interpo- 
lation of the EW(Pa-lines) using synthetic spectra based on 
the unblended Pal2, Pal4, and Pal7 lines (Fig.|6| We used 
standard LTE Castelli-Kurucz models (Castelli & Kurucz 
2008) to determine fluxes at different wavelengths corre- 
sponding to lines chosen as accretion indicators for all pro- 
gram stars. The values of their stellar masses A/* and radii 
i?* were taken from Table [2]to calculate luminosities in the 
line-indicators. 

The procedure of measuring the emission Balmer jump 
ADb is illustrated in Fig. fusing HD 100546 as an ex- 
ample. The ratio r of the initial unreduced spectra of the 
target (HD 100546) and the standard star of spectral type 
AO (HD 100928) is approximated polynomially between 
4000 A and 4600 A. The polynomial is extrapolated to 
3640 A and the ratio of ordinates r2 and r-y allows to de- 
termine the value of ADb = 2.5 log (r^rx). This procedure 
is identical to that used in the work of Donehew & Brittain 
d20TTT) . 



4 Spectral accretion indicators 

We examined 13 spectral diagnostics of M acc that have pre- 
viously been used in studies of TTS and BD. The basic in- 
dicator was the emission Balmer jump ADb, calculated by 
the method described in Sect. [3] We used the relation be- 
tween ADb and M acc presented in Fig. Q] Eight indica- 
tors, based on the measurements of I/ acc include Ha, H/3, 
He I A5876, Br 7 , Na I D, Pa/3, Pa 7 , and Ca II A8542. We 
considered also three indicators where flux in an emission 
line was used for the direct determination of M acc accord- 
ing to eq. 3: He I A5876, Ca II A8542, and Ca II A8662. 
Further in the text we use the abbreviations "(E)" and "(F)" 
to mark indicators based on the measurements of L acc and 
flux, respectively (e.g., Pa 7 (L), Ca II A8662 (F), etc.). The 




3600 3800 4000 4200 4400 4600 

A(i) 

Fig. 3 Spectral energy distribution (SED) in the ob- 
served spectrum of HD 100546 and of the normal main- 
sequence star HD 100928 of spectral type A0V in the 
range of AA3600-4600A covering the Balmer jump 
(top panel). The bottom panel illustrates the ratio r = 
SED(HD 100546)/SED(HD 100928). The solid line marks 
the polynomial approximation of r(A) between A4000A 
and A4600 A, and the dotted line shows its extrapolation up 
to A3640 A. The formula for the calculation of the emission 
Balmer jump is also inserted. 

last indicator is Hal0%, which was calibrated by Natta et 
al. (120041) . All these diagnostics are presented in Table 4. 

Selected regions of the observed spectra of the program 
stars are presented in Figs. |4]to|7] 



4.1 Ha emission 

The emission Ha profiles vary strongly from object to 
object (Fig. . HD 150193 and HD 190073 demonstrate 
PCyg-type profiles indicating the existence of a stellar 
wind between the star and the observer. The Ha profiles of 
HD 100546 and HD 101412 are double-peaked with a cen- 
tral absorption shifted relative to the stellar rest wavelength 
of the line. The Ha profile of HD 97048 is double-peaked 
too, but of a complex structure, and single emission profiles 
are observed in the spectra of HD 135433B and PDS 2. 
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Fig. 4 Normalized Ha profiles in the spectra of studied HAes. For HD 135344B, HD 150193, HD 176386, HD 190073, 
and PDS2 the displayed profiles are observed on dates MJD55357.329.371, 55352.047, 55465.138, 55446.078, and 
55410.343, respectively. Profiles of HD 135344B and PDS2 are shown in comparison with synthetic profiles calculated 
with the computer code SYNTH+ROTATE of (Pi skunov [19921 . The LTE models of Castelli & Kurucz d2"008"b are used for 
the stellar parameters listed in Table [2] In the case of PDS 2, the value v sin i = 30km/s was used to calculate the model 
profile. 
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Table 4 Mass accretion rate diagnostics that have been used in previous studies of TTs and BDs and are analyzed in this 
work 

Reference Diagnostic Mass range [M,/Mq] 

Fang et al. (2009) Ha (L), H/3 (L), He I A5876 (L) 0.1-2.0 
Herczeg & Hillenbrand ( [20081 Na I D (L), Ca II A8542 (F), Ca II A8662 (F), He I A5876 (F) 0. 1 - 1.0 
Dahm i2008l Pa/? (L), Pa 7 (L), and Ca II A8542 (L) 0.5-2.0 
Natta et al. <2004l Pa/3 (L), Hal0% 0.01 - 0. 1 
Gatti et al. (120081 Pa/3 (L), Pa7 (L) 0.1 - 1.0 



As can be seen in Fig. |4] the Ha emission is practi- 
cally invisible in the spectrum of the B9e star HD 176386. 
This fact was already discussed by other authors. Bibo et 
al. ([1992 ) characterized this object as a higher-mass analog 
of so-called weak-line TTauri stars (WTTS) with already 
dispersed accretion disks. On the other hand, Grady et al. 
( 1993 ) report that signatures of a matter infall onto the star 
is observed in the UV spectrum of HD 176386. If the disk 
of HD 176386 is already dispersed, calibrations derived for 
the magnetospheric accretion scenario are no longer valid 
because the balance of the emission from different parts of 
the circumstellar components (magnetosphere, disk, wind) 
changes. In any case, the standard model of magnetospheric 
accretion from the disk can hardly be applied to this object. 
Therefore, any empirical spectral calibration based on the 
MA model cannot be used. 

4.2 Other emission lines 

Fig. [5] illustrates the spectra of the targets in the region of 
the He I A5876 and Na I D lines, and Fig.|6]in the region of 
the IR Ca II triplet. The emission profiles of He I A5876 are 
in general single-peak or double-peak type in the spectra of 
the HAeBes, but in some cases a deep redshifted absorption 
overlaps the red emission wing (see for more details Bohm 
& Catala ll 9951 Mendigutia et al. 1201 lbt . An example of 
such a feature is seen in the spectrum of HD 97048 (Pig. O. 
This line originates in the highest temperature regions of the 
CS envelope involving the inner boundary of the accretion 
disk interacting with the magnetosphere, the matter streams 
infalling onto the star, and the region of the impact of the 
streams on the stellar surface. These regions are rather com- 
pact in size and screening of the stellar limb by the accreted 
flows can be an important factor in forming the observed 
line profile. This screening leads to the appearance of a red- 
shifted absorption component of the profile. Therefore, we 
have to be careful in using the He I emission for the M acc 
determination. At the time when the infalling stream screens 
the stellar limb, the M acc measurement can be considerably 
underestimated. 

The types of the Na I D line profiles (Fig. [5]l as well 
as of the IR Ca II triplet profiles (Fig. |6]l in the spectra of 
the targets are very diverse. Even in objects of similar spec- 
tral types (A2-B9) are the intensities of these emission lines 
completely different. The emission in the Na I D lines is 
clearly visible and the EW can be measured in three objects 
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Fig. 6 Normalized spectra of HD 97048, HD 100546, 
HD 150193, and HD 190073 in the spectral region contain- 
ing near-IR Ca triple lines. The positions of the Pa lines 
are indicated in the upper panel, and the positions of the 
Ca II lines in the lower panel. Strong emission near-IR Ca 
triplet lines appear only in two HAe stars, HD 150193 and 
HD 190073. 
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Fig. 5 The same as in Fig. |4] but for the spectral region around He I 5876 and Na I D. 




5900 



(HD 190073, HD 150193, and HD 101412). As for the Ca II 
triplet, the EW can be measured only in HD 190073 and 
HD 150193. It is remarkable that in our sample HD 190073 
is the star with the strongest intensities of the Na I D and 
IR Ca II lines. This object demonstrates the clearest devel- 
oped P Cyg-stmcture of the Ha profile (Fig.|4]l. The second 
object with an Ha profile of PCyg-type is HD 150193. It 
shows also rather intense emission in Na I D and IR Ca II. 
This fact leads us to the hypothesis that the appearance of 
considerable emission in the Na I D and IR Ca II lines is re- 



lated to the presence of a wind zone between the star and the 
observer. The signatures of the matter outflow are clearly 
seen in the Na I D lines of HD 190073 and HD 150193 in 
form of the P Cyg-structure. 

4.3 Balmer jump 

Fig. [7] illustrates the normalized spectra of the targets in 
the region of the Balmer jump in comparison with spec- 
tra of non-accreting stars of similar spectral types. Castelli- 
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Fig. 7 Normalized spectra of our targets in the region near 3600-4400 A. The normalization was made to the flux F(X) 
at the level of the model continuum at A = 4000 A. Spectra of standards of similar spectral type are given for comparison 
(dashed lines). 



Kurucz LTE (Castelli & Kurucz 2008 ) models were applied 
for normalization of the observed spectra of non-accreting 
stars to flux F\ at 4000 A at the continuum level. The r - 
dependencies (see Sect.[3]and Fig.0 were used to construct 
the normalized spectra of the targets. 

The EWs of all studied lines, the Hal0% width and 
ADb are listed in Table [5] together with the observing 



dates and magnetic/rotation phases calculated according to 
the ephemerides presented in Hubrig et al. (2011a 2 01 lbt . 
Mean accuracies of the measurements are presented in Ta- 
ble[6j The errors of the individual estimates of the equivalent 
widths were computed using a program that takes into ac- 
count the width and the maximum intensity of an emission 
line as well as the S/N ratio of the spectrum near the mea- 
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Table 5 Observing dates, magnetic/rotation phases calculated according to the ephemerides presented in Hubrig et al. 
(201 la, 201 lb), and measurement results of EWs, M / (HalO%), and ADb for HAes in our sample. 



Uuject 


MJL) 


Magn. 
phase 


T3 mi 


Do R 
rSLfj 


Pa^y 


C TT 

ca ii 

ODOZ 


Cm r £ n 
bw [A J 

/"., TT 

ca ii 


l in 


Mo t n 


IT, I 

ne i 

JO /D 


tip 


XX/ / TT _ 1 f\Of \ 

LKnvsj 




HD 97048 


55352.000 


0.43 


10.30 


21.90 


9.12 


- 


- 


30.8 


- 


0.0525 


1.58 


604 


0.25 


HD 100546 


55352.022 




16.05 


26.41 


9.72 


- 


- 


30.7 


0.035 


0.162 


2.55 


540 


0.18 


HD 101412 


55351.992 


0.71 


5.79 


7.77 


3.63 


- 


- 


13.6 


0.340 


0.145 


1.69 


560 


0.225 


HD 135344B 


55329.362 




2.27 


3.88 


1.47 


- 


- 


10.3 


- 


0.43 


1.57 


496 


0.15 




55352.034 




3.21 


4.10 


1.72 


- 


- 


14.5 


- 


0.32 


2.16 


501 


0.15 




55357.251 




2.65 


3.21 


1.39 


- 


- 


11.1 


- 


0.26 


2.00 


501 


0.15 


HD 150193 


55329.371 


0.59 


8.83 


7.66 


4.93 


0.87 


1.61 


11.9 


0.253 


0.530 


1.15 


517 


0.19 




55352.047 


0.81 


10.19 


9.80 


6.06 


1.06 


2.13 


13.6 


0.479 


0.479 


1.50 


561 


0.19 


HD 176386 


55329.387 


0.09 


0.392 


0.188 


- 


- 


- 


0.167 


- 


- 


- 


- 


- 




55463.171 


0.88 


0.184 


0.144 


- 


- 


- 


0.159 


- 


- 


- 


- 


- 




55465.138 


0.06 


0.270 


0.203 




- 


- 


0.175 


- 




- 


- 


- 




55466.166 


0.21 


0.231 


0.153 


- 


- 


- 


0.161 


- 


- 


- 


- 


- 


HD 190073 


55329.397 




9.97 


15.8 


8.04 


5.50 


6.58 


23.9 


1.94 


0.213 


2.64 


- 


0.19 




55410.250 




10.67 


18.8 


9.45 


5.91 


7.29 


26.8 


2.25 


0.337 


3.54 


- 


0.27 




55446.078 




10.97 


19.1 


9.74 


5.48 


6.87 


25.2 


2.14 


0.471 


3.58 


- 


0.27 




55463.182 




8.91 


18.1 


8.36 


4.65 


5.93 


20.5 


1.83 


0.364 


2.37 


- 


0.185 




55465.159 




8.65 


17.2 


8.67 


4.73 


5.92 


20.5 


1.82 


0.257 


2.56 




0.19 




55468.141 




8.40 


16.4 


7.82 


4.53 


5.56 


20.5 


1.79 


0.317 


2.18 




0.19 


PDS2 


55375.418 




1.60 


0.627 


0.358 






3.16 




0.0286 


0.34 


380 






55408.317 














5.69 




0.0745 


0.68 


348 






55409.410 














5.52 




0.0585 


0.54 


380 






55410.343 














4.50 




0.0464 


0.32 


396 






55411.394 




2.85 


2.05 


0.829 






5.90 




0.1068 


0.74 


365 






55446.307 




2.33 


1.00 


0.450 






3.30 




0.0621 


0.28 


424 






55463.210 




2.77 


1.30 


0.627 






5.25 




0.0734 


0.55 


425 






55464.208 




1.79 


0.95 


0.460 






6.10 




0.0335 


0.33 


409 





Table 6 Mean errors of measurements presented in Table |5] 



Object 


Br7 


Pa,S 


Pa7 


Ca II 

8662 


±cr(EW) [A] 

Call 
8542 


Ha 


NalD 


He I 

5876 


II 


W(Hal0%) 
[km/s] 


AD B 
[mag] 


HD 97048 


0.25 


0.40 


0.25 






0.45 




0.015 


0.10 


20 


0.03 


HD 100546 


0.30 


0.40 


0.25 






0.45 


0.015 


0.020 


0.15 


20 


0.03 


HD 101412 


0.20 


0.25 


0.20 






0.30 


0.05 


0.030 


0.10 


20 


0.03 


HD 135344B 


0.15 


0.15 


0.10 






0.25 




0.045 


0.10 


20 


0.03 


HD 150193 


0.25 


0.25 


0.20 


0.07 


0.10 


0.25 


0.04 


0.060 


0.10 


20 


0.03 


HD 176386 


0.08 


0.05 








0.05 












HD 190073 


0.25 


0.30 


0.25 


0.15 


0.20 


0.40 


0.10 


0.04 


0.15 




0.03 


PDS2 


0.10 


0.10 


0.06 






0.20 




0.02 


0.06 


20 





sured line. Further, the program considers the contribution 
of the photospheric background and the telluric line subtrac- 
tion. Since we used for the determination of ADb a method 
similar to that used by Donehew & Brittain (201 li, our er- 
ror estimates are of the same order as presented by these 
authors, about 0.03 mag. The errors for the individual mea- 
surements differ from the mean values presented in Table |6] 
by not more than 20%. 



5 Accretion rates 

5.1 Testing existing empirical calibrations 

We computed the M acc values for HAes in our sample using 
the measured quantities from Table |4] for indicators speci- 
fied at the beginning of Sect. [4] 

HD 176386 was excluded from our analysis for the rea- 
son discussed in Sect. 14.11 The formal calculation of M acc 
of this object leads to log M acc = —8.70 for B1-7, —9.03 
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for Br/3, and —9.63 for Ha. Other indicators are absent in 
all observed spectra of the star or too weak to be measured 
precisely. 

The indicator AD b was used only for six objects with 
-Mace > 10 _ Mq/jt . As was concluded in Donehew & 
Brittain (1201 j) , this diagnostics is inefficient for ADb < 
0.1 m . The width AV r (HalO%) cannot be measured cor- 
rectly for HD 190073 because a strong P Cyg-type structure 
overlaps the blue wing of the emission Ha profile. Taking 
into account that such features are not common to all HAes, 
there is a doubt if it is reasonable to use them at all. Still, it 
is of interest to compare the M acc estimations derived from 
these diagnostics with values obtained from other indica- 
tors. 

5.1.1 Criterion of applicability 

Our aim was to choose diagnostics that lead to results con- 
sistent with each other, where the basic indicator ADb is 
the only tracer connected with Af acc directly on the basis of 
model calculations. We introduce the criterion of applica- 
bility of a diagnostic "X" for measuring the M acc of HAes 
as follows: the indicator is applicable if the mean value of 
< log M acc (AD B ) - log M acc (X) > (further in the text as 
" ADb — X") is consistent with zero within the 1 er m errors, 
calculated for our sample. Here a m =a/y/n — 1 is the stan- 
dard error of the mean value and n is the number of objects. 

As a result, we identified a number of diagnostics sat- 
isfying this criterion (see Table [7). A small disagreement is 
present in the Pa/3 (Dahm l2008l) indicator: +0.30±0.09dex, 
and in the Ca II A8542(L) (Dahm 120081 indicator. In any 
case, the mean spread of all these differences type "ADb - 
X" is of the order of ±0.10-0. 15 dex. It is remarkable that 
if the lines Na I D and IR Ca II A8542 are clearly visible in 
the spectrum of an object, the M acc estimates derived from 
these indicators are in a satisfactory agreement with values 
obtained from other spectral diagnostics. 

The three remaining indicators of M acc that make use 
of F-type calibrations (Ca II A8542, Ca II A8662, and He I 
A5876) were additionally examined because the "ADb — 
X" values of two Ca II (F) indicators demonstrate differ- 
ences at larger amplitudes than all other tracers presented 
in Table |7j (+0.44 ± 0.32 dex and -0.28 ± 0.35 dex for 
Ca II A8542 (F) and Ca II A8662 (F), respectively). Further, 
the He I A5876 tracer shows a notable systematic shift of 
"AD B - X"= +0.69 ± 0.19dex. 

5.1.2 Testing F-calibrations 

The analysis of Af acc estimates that were obtained using 
these diagnostics demonstrate considerable systematic in- 
consistencies between each other as well as between F- and 
L-type calibrations for the same lines. These discrepancies 
are illustrated in Fig. [8] 

The estimates derived from Ca II A8542(F) are sys- 
tematically 0.75+0.03 dex higher than those from Ca II 
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Fig. 9 Comparison of the results of the log M acc [dex] de- 
termination for our targets obtained from the different Pa (3 
and Pa7 calibrations defined in Dahm (2008) and Fang et al. 
(2009). The dashed line indicates the equality of measured 
values. 

A8662(F) (top left panel). These two lines originate in the 
same circumstellar region, and we conclude that conditions 
in this region in HAes are different from those in TTS for 
which these two F-calibrations have been introduced. 

The mean difference between values obtained from 
A8542(F) and A8542(L) is +0.90dex for HD 190073 and 
+0.40 dex for HD 150193 (top right panel). A similar pic- 
ture is observed if we compare the M acc derived with He I 
A5876(F) and He I A5876(L) calibrations (bottom panel). 
The values obtained with the F-type calibrations are sys- 
tematically higher with the spread lying between 0.15 and 
0.70 dex for different objects. Each star has its own offset 
which is the same at all observed epochs pointing at some 
systematic factor. Based upon the results of this test, we as- 
sume that the L-type calibrations are best suitable for the 
HAes in our sample. We suggest that the observed system- 
atic differences between the values derived with the F- and 
L-type calibrations are likely to be caused by an insufficient 
accuracy of the c and d coefficients in the F-type calibra- 
tions applied to HAes, as well as by uncertainties in M* or 
i?* of the targets, used for the calculation of L acc . 

5.1.3 Other factors leading to systematic effects 

The uncertainties in stellar parameters (especially in the 
stellar radius i?*) can be a source of considerable system- 
atic errors in the M acc determination. As an example, we 
computed the mass accretion rate of HD 190073 separately 
with two different values of its stellar radius: R*/Rq = 3.6 
(Catala et al. 2007) and the more recent estimate R^/Rq — 
2.1 (Montesinos et al.|2009). In the first case, we have ob- 
tained a mean value of the accretion rate —6.23, with a 
spread at the level of standard deviation of ±0.17 dex us- 
ing eight spectral indicators derived from luminosities of the 
emission lines (L-type). This value strongly deviates from 
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Table 7 List of diagnostics satisfying our criterion of applicability (besides of ADb)- The error of the mean value a„ 
indicates the accuracy of determination. 



Diagnostic 

"X" 



Reference 



Mean value of 
"ADb - X" 



Number of 
objects 



[dex] 



Br 7 (L) 
Pa/3(L) 
Pa 7 (L) 
Ha(L) 
H/3(L) 

He I A5876(L) 
Na I D(L) 
Ca II A8542(L) 
Ha 10% 



Donehew & Brittain ( 12011b 

Dahm i2008l 

Dahm i2008l 

Dahm l2008"t 

DahmUQOD 

Fang et al. d2009t 

Herczeg & Hillenbrand ( [20081 

Dahm f2008l 

White & Basri 12003 1 



+0.19 ±0.15 
+0.30 ± 0.09 
+0.11 + 0.15 
-0.02 ± 0.12 
-0.13 + 0.09 
+0.17 + 0.17 
+0.10 + 0.16 
-0.22 ± 0.07 
+0.08 + 0.10 



■7.0 h 
-7.2 



C\2 

to 

to -7.4 

CD 



■7.6 h 
-7.8 



TV 



HD190073 



HD150193 



■6.4 - 
■6.6 



C\2 

5 -6.8 

CO 



-7.0 h 

-7.2 



n 1 1 1 r 

HD190073 



O HD150193 



-7.2 -6.9 -6.6 -6.3 
Call 8542 (F) 



-7.6 -7.4 -7.2 
Call 8542 (L) 



CO 

r- 
co 

LO 

I— I 

CD 

EC 



-6.0 r 

-6.5 r 

-7.0 r 

-7.5 r 

-8.0 r 

-8.5 ' r „ 
-9.0 



-I 1 1 1 - 



-i — I — I — r- 



I 1 1 1 1 I 



9 



J_ 



_1 I I L_ 



_1 I I L 



-9.0 -8.5 -8.0 -7.5 -7.0 

Hel 5876 (L) 



-6.5 



-6.0 



Fig. 8 Illustration of systematic discrepancies between the results of the log M acc [dex] determination obtained from 
comparing three pairs of diagnostics: Ca II A8542 (F) vs Ca II A8662 (F) (left top panel); Ca II A8542 (L) vs Ca II A8542 (F) 
(right top panel); and He I A5876) (L) vs He I A5876) (F) (bottom panel). The dashed lines indicate the position of equality 
of abscissa and ordinate. Results for different objects are marked by different symbols. In the upper panels the filled 
circles mark the six measurements for HD 190073, while open circles refer to measurements obtained for HD 150193. 
Ca line diagnostics can be used only for these two stars. In the bottom panel HD 97048 is marked by the open circle, 
HD 100546 by the open triangle pointing downwards, HD 101412 by the open triangle pointing upwards, HD 135344B by 
filled triangles pointing upwards, HD 150193 by open squares, HD 190073 by filled circles, and PDS2 by filled triangles 
pointing downwards. 
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the —7.14 obtained from ADb which, according to Fig.Q] 
is independent of M* and i?* of a star. Using i?* = 2.1 Rq 
leads to the mean value — 7.15±0.17, that is in good agree- 
ment with the estimation from ADb- 

Applying different empiric relations between Pa/3 and 
Pa7 luminosity, respectively, and the mass accretion rates 
that have been derived for the TTS by various authors to 
our sample, we find considerable and systematic discrep- 
ancies. This is demonstrated in Fig. [9] where we show the 
Af acc values for the targets that result from the calibrations 
by Gatti et al. d20Q8b versus those from Dahm (120081 1 . As 
can be seen in Fig. [9] the difference between the two cali- 
brations is practically absent for small log Af acc (less than 
—8.0) and becomes up to 0.40 dex for log M acc > —7. We 
have to mention that Gatti et al. ([2008 ) have calibrated their 
relation over the logA/ acc range from —9.6 to —8.2 (M* 
range 0. 1-1.0 Mq) and Dahm d2"008i l from -8.7 to -7.2 
(M, range 0.5-2.0 Mq), that is more similar to values M acc 
and M* for our targets. Therefore, we used in our study the 
diagnostics Pa/3(L) and Pa7(L) introduced by Dahm (2008 ). 

5.2 The deduced accretion rates for the targets in our 
sample 

Fig. [10] illustrates the mass accretion rates obtained for 
seven program stars (except HD 176386) that have been de- 
rived using the indicators that passed successfully our appli- 
cability test (see Sect. 15.1. lb . Because the F-type diagnostics 
lead to systematic differences in the M acc estimates in com- 
parison with the values obtained with the other indicators 
(see Sect. l5.1.2l and Fig. [8]), we did not include the values of 
M acc derived from F-type indicators in the calculations of 
the mean values and the uncertainties of the mass accretion 
rates. As a result, only L-type indicators were used (except 
Hal0% and ADb) and, further in the text, we dropped the 
abbreviation (L) in the names of the indicators. 

The mean values of M acc over the observed epochs and 
over several diagnostics as well as their standard devia- 
tions are presented in Table [8] where estimates from Gar- 
cia Lopez et al. (2006) and Mendigutia et al. ( 1201 1 al > are 
also given for comparison. One can see that, in spite of the 
fact that the observations of Garcia Lopez et al. (20061 were 
carried out in 2004, the discrepancies are within the limits 
of the standard deviations obtained in both works. On the 
other hand, the results obtained in Mendigutia et al. (201 la) 
demonstrate significant differences with our estimates. The 
difference is +0.94dex for HD 150193 and +2.15dex for 
HD 190073. Using data on A v , L(Ha), and L(Brj) pre- 
sented in Mendigutia et al. (1201 1 al > and from our paper, we 
estimated contributions of various factors forming the total 
discrepancy. These contributions are: +0.2 — 0.3 dex fol- 
low from differences in fluxes measured on different dates, 
+0.4 — 0.5 dex result from distinctions in A/ acc calibrations 
(for both HD 150193 and HD 190073). The contributions 
from differences in adopted stellar parameters M* and i?* 
are +0.2 dex and +1.5 dex for HD 150193 and HD 190073, 



respectively. Such a significant discrepancy for HD 190073 
is mainly related to the very large value of the stellar radius, 
R* = 8i?Q, used in Mendigutia et al. (1201 lab . This radius 
is by a factor of four larger than the value adopted in our 
work (Table |2). 

We note that the existence of a HAeBe with T e s 
less than 10000 K, Af* = 5.1 Mq, and i?„ = 8i? 
(Mendigutia et al. 1201 1 al > is not consistent with the PMS 
evolution model by Palla & Stahler (1993). According to 
their model, a PMS object with the initial Af* = 5 Mq 
starts to become visible at the birthline with T c g of about 
12 000 K, and later, closer to the end of the PMS stage near 
the main sequence, its effective temperature will become 
about 16 000 K. With the age of about 10 6 years given for 
HD 190073 in Hubrig et al. ( 120091 and Mendigutia et al. 
(2011a), the star has to be located in the H-R diagram al- 
ready near the main sequence, and with M* = 2 — 3 Mq 
(Table [5} its stellar radius i?„ cannot be as large as 8 Rq. 

As was discussed in Sect. 12.21 the M acc calibrations de- 
rived in Mendigutia et al. (1201 1 at lead to very large values 
of the mass accretion rates for a considerable part of their 
targets, which are not in agreement with the masses of ac- 
cretion disks (0.01 — 0.1 Mq) and the length in time of the 
PMS stage for HAes (10 5 — 10 6 years). They demonstrate 
also significant differences with the results of Donehew & 
Brittain (1201 11 1 and Garcia Lopez et al. (2006). 

Our M acc determination is based on: 

- consistence of values derived from several independent 
calibrations, 

- agreement of our estimates with those obtained in Gar- 
cia Lopez et al. (2006]!, and 

- accordance of our results with the predictions of the 
PMS evolution model by Palla & Stahler JT9931 ). 

Clearly, for HAEs with accretion disk masses of the order 
of 0.01 - 0. 1 Mq and a length of the PMS stage of HAEs of 
10 5 — 10 6 years, the expected mean accretion rates during 
the PMS stage are of the order of 10~ 6 - lO _8 Af /yr. 
These expected rates are similar to the values obtained in 
our work. 

Comparing the accuracy of our estimates with the re- 
sults obtained by other authors, we can see that the differ- 
ences are not significant (see Table |9). 

6 Variability of the mass accretion rates 

Four out of seven of our targets were observed on more than 
one occasion (from two to eight). As was discussed in the 
previous sections, all emission lines chosen as indicators of 
Af acc are formed not only in the magnetosphere, but also in 
the disk and in the wind, and their contribution to the en- 
tire emission is different for each line. Investigating a tem- 
poral behaviour of Af acc using different spectral diagnos- 
tics, we have to be certain that a variability revealed from 
these indicators relates to a change of Af acc and not to varia- 
tions of physical and kinematical conditions in the disk and 
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Fig. 10 Mass accretion rates of the targets determined from different indicators, as labeled. The solid line shows the 
mean log M a cc value and the dashed lines indicate the region of ±ler uncertainty where a is the standard deviation. 



Table 8 Mean mass accretion rates. The error is given as the standard deviation a. Literature data on accretion rates are 
taken from Garcia Lopez et al. (f2006) and Mendigutia et al. (1201 lal l. 



Object 


Number of 


Number of 


log Macc 


log M a cc 


log Macc 




indicators 


observations 


[M /yr] 


[M /yr] 


[M /yr] 










Garcia Lopez et al. (120061 


Mendigutia et al. d201 lab 


HD 97048 


8 


1 


-6.80 ±0.38 


-7.17 




HD 100546 


8 


1 


-7.23 ±0.13 






HD 101412 


9 


1 


-7.04 ±0.15 






HD 135344B 


8 


3 


-7.69 ±0.30 


-8.27 




HD 150193 


10 


2 


-7.06 ±0.20 


-7.29 


-6.12 


HD 190073 


9 


6 


-7.15 ±0.16 




-5.00 


PDS2 


7 


8 


-8.68 ±0.36 
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Table 9 Mass accretion rate determinations using several diagnostics presented in this work and the results published in 
the literature. 









Mean standard 


Error of the 


Paper 


Range of log M acc 


Range of individual 


deviation a 


mean value a m 




[M /yr] 


deviations [dex] 


[dex] 


[dex] 


This paper 


-8.68 + -6.69 


0.13-0.36 


0.25 


0.11 


Dahm <|2008l 


-8.72 ^ -7.20 


0.1 - 1.4 


0.30 


0.12 


Rigliaco et al. 12011) 


-9.86 


0.4-1.4 


0.45 


0.14 



the wind. For example, the emission in such lines as Ha, 
H/3, and DNa I originates predominantly in the wind, while 
that in the He I line and near the Balmer jump in the high- 
temperature region close to the stellar surface (Pogodin et 
al. l20121 l. If the variations take place particularly in the outer 
circumstellar envelope, they are reflected differently in each 
of these lines. As a result, the amplitudes and even the char- 
acter of the measured M acc variability derived from differ- 
ent spectral calibrations can be different too. In such a situ- 
ation, the revealed changes of M acc must be considered as 
an artifact. 

6.1 Temporal behaviour of individual accretion 
diagnostics 

Fig. QT| illustrates variations of M acc of the two targets 
HD 190073 and PDS 2 with the largest number of observa- 
tions determined from all indicators. Since the amplitudes 
of the variability are smaller than the discrepancies between 
the values obtained from different diagnostics, we analyzed 
not the values themselves but the residuals relative to the 
mean value for each indicator. As one can see in Fig. QT| 
both objects demonstrate variability which: (a) correlates in 
a majority of the indicators, and (b) shows amplitudes of 
the variations that are very similar for different diagnostics. 
It could be intrinsic M acc variability but also some geomet- 
ric effect that affects in the same way our view of all the 
different emitting regions. 

However, one indicator exists for both objects showing 
a distinct temporal behaviour of M acc : it is He I A5876 for 
HD 190073 and Hal0% for PDS 2. The character of the 
variability of the Af acc value from the He I line in the spec- 
trum of HD 190073 is similar in comparison with other in- 
dicators, but shows a much larger amplitude. In the case 
of PDS 2, the estimations derived from the Hal0% width 
strongly differ from all others. A possible cause for this can 
be related to the existence of small-amplitude variations tak- 
ing place in the circumstellar gas outside the magnetosphere 
(the outer disk and the wind) which may influence the full 
width of the profile but might be insufficient to have a sig- 
nificant effect on the EW of Ha. This could explain a rather 
small distorting influence of the non-magnetospheric vari- 
ability onto the M acc derived from the Ha(L) diagnostic and 
a rather significant - onto that obtained with the Hal0% in- 
dicator. 



The other two stars with more than one available spec- 
trum, HD 150193 and HD 135344B, demonstrate the tem- 
poral behaviour of M acc derived from all calibrations that is 
similar to the case of HD 190073. The only exception is the 
He I A5876 diagnostic. Variability of Af acc values derived 
from this indicator is quite different from that obtained from 
all other indicators. Presently, no strict model for the forma- 
tion of the He I A5876 line in the high-temperature circum- 
stellar gas around HAeBes exists. We can only assume that 
a multi-component variability is likely taking place that is 
connected with a change of geometric configuration, optical 
thickness and emissivity of the gaseous flows accreted onto 
the star inside the magnetosphere. But this has a rather small 
amplitude and the mean value derived from He I differs not 
so strongly from all the others estimates (see Fig.fToV 

6.2 Global accretion variations 

As a measure for the overall accretion variability, Fig. [T2lil- 
lustrates the temporal behaviour of the mean values of M acc 
for the four targets derived from several diagnostics. The er- 
ror of the mean value a m =al\/n — 1 (where a is the stan- 
dard deviation and n is the number of the used diagnos- 
tics) is shown by vertical bars. Indicators showing a tem- 
poral behaviour of M acc different from all others were ex- 
cluded from the calculation of the mean values and their 
errors (these are He I A5876 for HD 190073, HD 135344B, 
and HD 150193, and HalO% for PDS 2). 

HD 190073 demonstrates variability on the timescale of 
tens of days with a spread of the values of about 0. 1 dex with 
an accuracy at la m level of ±0.01 dex. Short-term varia- 
tions (from night to night) are undoubtedly found in PDS 2. 
The observations of HD 135344B and HD 150193 were car- 
ried out only on two to three different nights, separated by 
several to tens of days. Thus, it is impossible to examine the 
cause of the variability. 

An effective way to analyze the amplitudes and the 
timescale of accretion variability was suggested by Nguyen 
et al. ( 120091 1. where a change in amplitude versus the time 
interval between observations was constructed for 40 clas- 
sical TTS. The authors sampled the timescale range from 
hours to months and concluded the maximum extent of the 
M acc variability is reached after a few days and amplitudes 
of variations are as a rule not more than 0.5 dex. This result 
has been confirmed by other authors. Costigan et al. (12012) 
bolster it with a sample of 10 accreting TTS observed in 



© 2010 WILEY- VCH Verlag GmbH & Co. KGaA, Weinheim 



www.an-journal.org 



Astron. Nadir. / AN (2010) 



17 



<D 



M 
O 



Hel 



I 1 1 1 1 I 1 1 1 1 I 

HD190073 

residuals 




□ 



u 

6 



I 0.10 



dex 



55300 55350 55400 



-i — i — i — r- 



□ 

Hal0% 



PDS2 
residuals 




u 

d 
o 

■•d 



□ 



I 0.20 dex 



J_ 



55450 55350 
MJD 



55400 



55450 



Fig. 1 1 Temporal behaviour of log M acc in HD 1 90073 and PDS 2 obtained from different indicators, constructed relative 
to the mean value for each indicator. Thin solid lines join neighbour points for each indicator for better illustration. The 
temporal behaviour of the residuals log M acc obtained from the indicator He I (for HD 190073, left panel) and from the 
indicator Hal0% (for PDS 2, right panel) are marked by open symbols and thin dashed lines. Dotted lines indicate positions 
of zero for each indicator. 



the Longterm Accreting Monitoring Project (LAMP) show- 
ing that accretion variation amplitudes do not increase when 
extending the timescales of the monitoring beyond a week. 

In Fig. [T3] we present our analysis of the M acc vari- 
ability as a function of timescale following Nguyen et al. 
(2009). Our results are fully in agreement with those ob- 
tained previously by other authors. In the case of PDS 2, 
variations on a timescale of about tens of days with ampli- 
tudes of up to 0.4 dex are observed. Besides of that, variabil- 
ity on a timescale near one day is clearly seen with ampli- 
tudes of up to 0.3 dex. Such a timescale is comparable with 
magnetic/rotation periods expected for our targets (Hubrig 
et al. 120091 1201 lab . Therefore variability on such a short 
timescale is likely to be connected with a modulation of 
profiles and intensities of the spectral lines by a rotating 
magnetosphere. Using data from Table 2 (i?* and vsini), 
we estimated the rotation period of PDS 2 P = 2.7 sini days. 
This object is likely to be oriented relative to an observer 
close to "pole-on", and its inclination angle i is expected 
to be of order of 15°. Such orientation is confirmed by a 
single peak emission Ha profile (Fig. |4j (see Grinin & Ros- 
topchina [T9"9"6i l. 

Due to the insufficient number of observations, it is 
not possible to analyze the M acc variability on different 
timescales for the other three targets. The only conclu- 



sion can be made that a variability has been revealed on 
timescales from several days to tens of days. In the case of 
HD 190073 we can derive the rotation period P = 9.1 sini 
days using data from Table 2. In Hubrig et al. (2009]! the 
upper limit of i is given as i < 40°. A PCyg-type of the 
Ha profile (Fig. [3]l evidences that this object has to be of 
an intermediate orientation relative to an observer (Grinin 

6 Rostopchina 1 1 996t . With the expected value of i ~ 40° 
we estimate P^6 days. Since no observations separated by 
short time intervals have been carried out for this target so 
far, we cannot examine a presence of rotational modulations 
of the line profiles in the spectrum of HD 190073. 

Another plausible alternative could be that the variabil- 
ity is the result of stochastic variations in the accretion 
rate through the disk, as one might expect if the magneto- 
rotational instability is the primary source of viscosity in 
accretion disks. More spectral data obtained on different 
timescales (from days to tens of days) are needed to rig- 
orously test the character of the variability of the Herbig 
Ae/Be stars. 

7 Conclusions 

We present results of observations of a sample of eight mag- 
netic Herbig Ae stars obtained with the X-shooter spectro- 
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graphs installed on the 8 m UT2 at the VLT. This spec- 
trograph covers simultaneously the whole spectral range 
from the near-UV to the near-IR (300-2500 nm), provid- 
ing medium resolution, high-quality spectra. We examined 
13 different spectral accretion diagnostics derived earlier on 
the basis of TTS and BD observations on their applicabil- 
ity to HAes. This applicability has been confirmed for a 
number of the indicators. The criterion of the applicabil- 
ity was based on a consistence of the results of Af acc de- 
termination one with another and with the only direct ac- 
cretion indicator ADb, based on the CS emission near the 
Balmer jump. We have used the diagnostics satisfying the 
criterion to compute Af acc of all targets, which are ranged 
from 2 x 10~ 9 to 2 x 10~ 7 Mq /yr with standard devia- 
tions of 0.25 dex relative to the mean value. The targets in 
our sample demonstrate a large variety of spectra. One of 
the eight targets, HD 176386, was excluded from our anal- 
ysis as an object that has a "non-developed" accretion disk. 
It was not possible to use all 13 indicators in all objects. 
The number of the indicators used varies from 7 to 10, de- 
pending on the target. An important factor limiting the ap- 
plication of all indicators is the necessity to re-calibrate a 
number of empiric relations constructed earlier for objects 
showing small Af acc < 10 _9 AfQ/yr. Their application for 
HAes with higher (M acc ~ lO _7 Af0/yr) can lead to sig- 
nificant systematic errors. We would like to emphasize that 
additional future spectral observations of a larger number of 
HAes are needed to improve all calibrations and to expand 



them to the sample of earlier-type PMS objects with higher 

Afacc- 

An important cause of errors is the uncertainty in deter- 
minations of A/* and especially f?* for a number of targets 
when the indicators are based on luminosity determination. 
For objects with small accretion rate values (log Af acc < 
—8.5) errors of measurements of EWs and widths of emis- 
sion profiles become an additional factor for the measure- 
ment accuracy. 

Our investigation of variability of Af acc has shown that 
all four Herbig Ae stars observed on more than one ocassion 
demonstrate a change of Af acc of about 0.10 -j- 0.40 dex on 
the timescale of tens of days. One object observed on a short 
timescale, PDS 2, shows also night-to-night variations with 
an amplitude of up to 0.30 dex. This variability might be re- 
lated to the modulation of spectral parameters by a rotating 
magnetosphere. In the case of even longer-term variability 
(r R3 tens of days) it remains unclear whether this variabil- 
ity is connected with a change in the accretion regime or is 
a result of rotational modulation. We note that due to the 
small number of observations of targets with determined 
magnetic/rotation periods, not much can currently be con- 
cluded on the role of magnetic fields in the dynamics of 
the accretion processes. Future multi-epoch observations of 
magnetic HAes with X-shooter will be extremely useful to 
better understand the nature and variability of Af acc spectral 
diagnostics to constrain more realistically magnetospheric 
accretion models for Herbig Ae stars. 
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Fig. 13 Amplitude of variations of the mean log M acc of 
four targets as a function of time elapsed between the ob- 
servations. Upper panel: Results forHD 190073 (filled sym- 
bols) and HD 135344B (open symbols). Bottom panel: Re- 
sults for PDS 2 (filled symbols) and HD 150193 (open sym- 
bols). 
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